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I review present challenges that QCD in extreme environments presents to lattice gauge theory. Recent data 
and impressions from RHIC are emphasized. Physical pictures of heavy ion wavefunctions, collisions and the 
generation of the Quark Gluon Plasma are discussed, with an eye toward engaging the lattice and its numerical 
methods in more interaction with the experimental and phenomenological developments. Controversial, but 
stimulating scenarios which can be confirmed or dismissed by lattice methods are covered. In the second half of 
the talk, several promising developments presented at the conference Lattice 2002 are reviewed. 



1. Introduction 

We live in interesting times. 

In the world of high energy and nuclear physics 
the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory is beginning its 
search into the new realm of high temperatures 
and low but nonzero chemical potentials. 

These experiments will surprise us. Experi- 
ments have a habit of doing that. They humble 
us. They will show us new directions. They will 
make this talk obsolete. 

I want to emphasize that Lattice Gauge The- 
ory, which has become an evolutionary rather 
than a revolutionary field, should take a more ac- 
tive part in these developments. It should rise to 
the challenge of real data to find new methods 
and ideas and extend its comfortable Euclidean 
ways to describe, predict and learn from scatter- 
ing experiments. Pioneering has the potential of 
reinvigorating the field. Lattice gauge theory has 
provided a solid estimate for the critical temper- 
ature to make the quark gluon plasma as well as 
not-so-solid estimates of the required energy den- 
sity, the magnitudes of screening lengths etc. But 
there is much more to predict. ..hopefully before 
the experiments... and doing so will be particu- 
larly exciting. And fun. 

I think that there arc promising areas for such 
developments and I will discuss some through- 
out this talk. Light-cone wavefunctions of nu- 
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clei, dense with gluons, quarks and anti-quarks, 
collisions with unexpectedly high multiplicities 
and signs of early development of a quark-gluon 
plasma will be reviewed. Along the way sug- 
gestions for additional or new lattice based re- 
search will be made. Other subfields of high en- 
ergy physics are already interacting with the data. 
Lattice gauge theorists should become more ac- 
tive participants. 

It hardly needs to be emphasized here that 
when Lattice Gauge Theory develops a method 
to analyze a problem, it can do so from first prin- 
ciples, with ultimately no approximations. Em- 
phasis on the word " ultimately" , because our lat- 
tices are yet too small, our quark masses are yet 
too large, our statistics are yet too few, etc. But 
the field is making steady progress on all these 
issues, helped particularly by its use of improved 
but practical Actions. Lattice gauge theory is the 
only approach that can handle the physics of the 
three crucial length scales of hadronic dynamics 
at one time, in a consistent calculation. At suf- 
ficiently short distances lattice calculations and 
simulations confirm asymptotic freedom, which is 
so essential to exposing the underlying gluon and 
quark degrees of freedom of QCD as experiments 
enter the quark-gluon plasma phase. At more 
moderate distances where the running coupling 
is in the intermediate range and semi-classical in- 
stanton configurations are breaking the anoma- 
lous U{1)a symmetry and, through properties of 
their ensemble, are breaking chiral symmetry and 
are producing the constituent quark masses, lat- 
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tice gauge theory is at its best elucidating the 
complex crossover physics of this range of length 
scales. Finally, at greater length scales, Lat- 
tice methods confirm confinement, the fact that 
hadronic states are color singlets and the confin- 
ing dynamics comes through thin, discernable but 
breakable flux tubes. QCD will not be fully un- 
derstood until these three qualitatively different 
ranges of phenomena are incorporated into one 
tractable analytic approach. Crucial hints needed 
to accomplish this will come from lattice studies. 

These are grand, and over-stated words. But 
progress is occurring. 

It is a pity that progress is not occurring on 
the challenge of producing a lattice simulation 
method for QCD at nonzero Baryon chemical po- 
tential fjL and vanishing temperature T . The in- 
famous sign problem of the fermion determinant 
continues to stand in our way. The same prob- 
lem has bedeviled condensed matter physicists for 
almost 40 years. 

A theme of this talk is that studying extreme 
environments teaches us how QCD works under 
ordinary conditions. I believe that when we have 
licked the sign problem, or have found a new for- 
mulation of nonperturbative QCD free of it, we 
will be at a new level in our understanding of 
how QCD really works and makes its low lying 
Baryonic excitations. 

In the second half of this review I will discuss 
several interesting contributions at this confer- 
ence. These will include progress in mapping out 
the low fi (chemical potential), high T (tempera- 
ture) part of the phase diagram of QCD, follow- 
ing the seminal work of Fodor and Katz. Real 
time spectral functions for the production of lep- 
ton pairs will be discussed as will the dispersion 
relation of pions below but near the transition 
to the plasma. A first step toward simulating 
phenomenologically interesting cutoff Four Fermi 
models of the transition will also be mentioned 
briefly. 

2. What's so Great about the Quark Gluon 
Plasma? 

As emphasized by E. Shuryak the mass 
scales of the Quark Gluon Plasma are different 



and, importantly, smaller than those of the more 
familiar hadronic phase. The hadronic phase 
breaks chiral symmetry, the quark gluon plasma 
does not. The hadronic phase confines quarks, 
the quark gluon plasma does not. The binding 
mechanism in the hadronic phase is nonperturba- 
tive while the screening mechanism in the quark 
gluon plasma is perturbative. 

We know from deep inelastic scattering that 
the substructure scale in the hadronic phase is 
Q-^ « 1 — 2 GeV. At this the running cou- 
pling is becoming large enough that perturbation 
theory is failing to give accurate estimates and 
nonperturbative effects are competitive. 

By contrast perturbative screening masses in 
the quark gluon plasma are Mf'"°" « 0.4 GeV. 
and Mt^"-^'' « 0.3 GeV. for temperature above 
but near T,. where the plasma first appears. 

The finer level spacings in the plasma act as 
a fine resolution grid to the dynamics in the 
hadronic phase. A collision which starts in the 
hadronic phase and ends in the plasma phase is 
especially sensitive to the levels and dynamics of 
the hadronic phase because of the relative wealth 
of open channels in the plasma phase. 

3. Overview of RHIC Collisions 

RHIC collisions can be analyzed in four stages 
which are thought to be characterized by distinct 
time scales. 

First, viewed from the center of momentum 
frame, there are right moving (positive pz) and 
left moving (negative Pz) nuclei. Each of them 
can be described by infinite momentum frame 
wavefunctions which give the probability ampli- 
tudes for each nuclei to consist of ensembles of 
quarks and gluons. The parton model, suitably 
improved with asymptotic freedom to describe its 
short distance features, should provide a frame- 
work for these initial states. 

Second, there is the collision between the con- 
stituents in each nuclei. Experiments will shed 
light on the typical momentum transfers and mul- 
tiplicities of these underlying collisions. 

Third, there is the development of a thermal- 
ized plasma of quarks and gluons. 

And fourth, there is the development of the fi- 
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nal states of hadronic debris from the hot soup of 
colliding and produced constituents. 

3.1. The Initial State 

Begin with the lightcone wavefunctions of the 
projectiles. The rightmover is, 

J Vi 
n i ' 

6^Y1 - E V^)\m,Kl; ?7„, i?„ > (1) 

Recall from infinite momentum frame |^ or 
light cone perturbation theory |Q that rji is the 
longitudinal fraction of the i*^ constituent. It is a 
positive variable, < ?7i < 1 and the constituents 
in each wavefunction ipn account for the longitu- 
dinal fraction of the projectile, r/i = 1. The 
absence of vacuum structure in Ivf > makes this 
wavefunction a useful formulation of the initial 
state of the collision. 

The bulk features of the wavefunctions of the 
projectiles are conveniently displayed on a rapid- 
ity plot, r = i ln[(po+Pz)/(po-Pz)]- The rapidity 
is particularly handy because it simply translates 
under a boost along the z axis. One can plot the 
density of partons of both projectiles on a rapid- 
ity plot and deep inelastic experiments at a fixed 
indicate that the plot is essentially flat and 
boost invariant and of extent Ins, where s is the 
usual scattering variable, four times the square of 
the center of momentum energy. 

Experiments indicate that parton rapidity dis- 
tributions have some simple, statistical features. 

1. Short range correlations in rapidity. The 
parton distributions lose memory, with a correla- 
tion length of 1 — 2 units, in rapidity. Screening 
in rapidity is an essential feature. Therefore, fea- 
tures like the quantum numbers of each projec- 
tile are limited to the edges of the rapidity plot. 
Unfortunately, this means that in very high en- 
ergy heavy ion collisions, most of the partons are 
not influenced by the Baryon number of the pro- 
jectiles. The collisions can be characterized by 
a high energy density and temperature but not 
high chemical potential fi. 

2. The "central region" and is universal and 
"vacuum like". The length of the central region 
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Figure 1. Rapidity distribution of partons in a 
nucleus-nucleus collision. 



grows like In s. 

Although first principle calculations of light 
cone wavefunctions are beyond us, models hav- 
ing these features have been known for decades. 
Multiperipheral graphs of QCD can populate the 
rapidity axis uniformly and have suggested that 
gluons are the most likely parton species in the 
central region ||^. 

Asymptotic freedom effects these considera- 
tions profoundly. Recall that in the naive par- 
ton model, partons interact only softly with one 
another, never exchanging high transverse mo- 
menta. But if we consider the structure of the 
wavefunctions at high spatial resolution where 
the running coupling is small and perturbative 
asymptotic freedom applies, then high transverse 
momentum exchanges and short distance addi- 
tional structure are inevitable Call the 
gluon distribution function in momentum space 
ryG(?7,Q^)- This is the distribution function that 
would be resolved by a probe with resolving power 
6Xj^ ^ -Qj. Models and theoretical analysis 
of scattering data suggests that rjG{irj,Q'^) ap- 
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proaches a constant for small and fixed . 
However, as is increased one discovers that 
gluons can split perturbatively into three or four 
gluons of lesser longitudinal fractions r]i. So, for 
sufficiently small 77, i]G{ri, Q^) must be an increas- 
ing function of . This is a general feature of the 
asymptotically free parton model Q and pertur- 
bative evolution equations make this simple phys- 
ical idea quantitative For quark distribution 
functions which are accessible to deep inelastic 
electron and neutrino scattering, these ideas are 
well tested phenomenologically. 

If we view the parton density on the plane 
transverse to the collision axis p^, then at low 

we see just a few gluons of size 5X\ ~ 
but as is taken larger we resolve many more 
gluons of smaller size. This progressive develop- 
ment has lead to the idea of "gluon saturation" 
in the light cone wavefunction of a nucleus . 

To understand this idea, consider the wave 
function and imagine that it is coming at you 
and you have a snapshot of the transverse plane. 
Call the transverse size of the nucleus Ra- Each 
parton of transverse momentum Q in the wave- 
function occupies an area tt/Q^, as suggested by 
the uncertainty principle. The parton's cross sec- 
tion is proportional to the square of the running 
coupling and its geometrical size, a ~ aslQ^)^- 
But partons will overlap in the transverse plane 
when their number Na becomes comparable to 
which is -^-jqjj- One supposes that when 
Na is this large, the overlap stops further growth 
in the density and " saturation" has been achieved 
1^. The transverse momentum of the partons at 
saturation is Ql ^ aJQ^)^ ~ A^^^. The num- 
ber of partons at the saturation level is Ns ^ 
^"(^Qi-^ ~ A and these are the number of par- 
tons of size scale 1/Qs which could material- 
ize in the quark gluon plasma. In addition 
sets the scale for the transverse momenta of par- 
tons in the central region of the rapidity plot, 

1 dN _ 1 p( K± \ 
a drd^K^ "sCQ^) ^ Qs 

These formulas can be made more precise by 
considering longitudinal fractions Q. The dis- 
cussion here is just a simplified, but, hopefully, 
intriguing introduction. When numbers for the 



various nuclei and the energies available at exper- 
imental facilities are substituted here, one finds 
that Q\ ^ \ — 2 GeV. This means that apply- 
ing perturbative QCD is suspect and competitive 
nonperturbative effects should be expected. This 
is potentially the realm of lattice gauge theory. 
The relationship of Euclidean lattice calculations 
and the Minkowski scattering formulation appro- 
priate to heavy ion collisions must be addressed. 

4. Space Time Evolution of the Collision 

Let's comment on the four different stages of a 
heavy ion collision separately. The physics issues 
in each are quite distinct. 

1. Wave functions approach on another. 

As reviewed above, the phase space parton den- 
sity is quite large ~ ^ {Q^) ' Classical and semi- 
classical estimates should apply. In particular, 
the color per volume is large, so the lack of com- 
mutivity, [Q", Q'^] = if^^^Q"^, appears to be a rel- 
atively small effect, and the gluon field can be 
treated classically for many purposes . 

2. Interactions. 

The gluon saturation picture suggests that the 
characteristic time of interactions is quite small, 
Unt ~ 0.1 — 0.3 fm./c and the energy 

density is very, very large, eint ~ a,{Q-^.) ~ 20 
GeV./fm.'^. Such a large energy density would 
place the system deep in the quark gluon plasma 
phase. This is good news! Estimates of tint in 
the past were much more modest, on the order of 
several GeV./fm.^, at most. 

3. Thermalization. 

After the collision the quark gluon matter ex- 
pands and thermalizes. The experimental data 
and models suggest that the time scale for this to 
occur is ttherm 0.5 — 1.0 fm/c. 

RHIC data, especially that showing the spatial 
distribution of particles produced in non-head-on 
collisions, strongly suggests that the interactions 
at early times are strong and the quark gluon 
plasma appears very early in the evolution of the 
final state 0. The basic parton-parton interac- 
tions must be strong and involve high multiplic- 
ities. Perturbative processes alone are not suffi- 
cient to explain the data. 
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4. Hydrodynamic expansion. 

The quark gluon plasma maintains thermaliza- 
tion and expands until decoupling sets in at tdec 
lOfm./c. In the expansion and production of the 
final state hadrons, the fastest particles are pro- 
duced last, as in the " inside-outside" cascade of 
electron-positron annihilation processes The 
expansion process is essentially one-dimensional, 
along the collision axis, until the latest stages of 
the development of the final state. A particle's 
rapidity turns out to be strongly correlated with 
the space-time rapidity, ^ ln[(a;o ^ x^^)/ (a;o ~ a;3)] > 
of its point of creation. 

5. Suggestions from RHIC and Evidence 
for a Quark Gluon Plasma 

RHIC experiments are hard and complicated 
and subject to many phenomenological models. 
Typically several models, each in stark contra- 
diction to the other, can describe a limited set of 
heavy ion data equally well. The "smoking gun" 
for the existence of the quark gluon plasma has 
yet to be found, but there are several pieces of 
data which are nicely explained, or at least, ac- 
commodated by it. Consider a few pieces of the 
puzzle. 

1. Jet Quenching. 

The heavy ion collisions show little or no sign 
of jets for Pi < 4 GeV. This is interesting in light 
of the fact that if you scale single particle hadron 
spectra from proton-proton data to A-A data, you 
overshoot the jet data very significantly Q. 

It has been argued that this is evidence for 
the existence of an extended space-time region 
of the quark gluon plasma. The idea is that 
two colliding gluons in the plasma annihilate 
into an energetic quark anti-quark pair with size- 
able tranverse momenta, but their interactions 
with the plasma medium saps the energy out of 
the energetic quark and its partner, eliminating 
("quenching") those jet-like characteristics. Of 
course, very high energy jets are uneffected by 
the medium because of the short range nature of 
the most prevalent interactions on the rapidity 
axis, and jets with energies much larger than 4 
GeV. are seen. 

2. Suggestions of Dcconfincmcnt and Chiral 



Symmetry Restoration. 

Lepton pairs have long been cited as effective 
probes into the internal dynamics of the quark 
gluon plasma. In particular, if the p and 4> existed 
in the plasma, then the leptonic decays, p — > ee 
and — > ee would be easily seen, as in proton-A 
collisions. These peaks are missing, however 0. 
It is tempting to interpret this experimental result 
as evidence that the light hadrons have "melted" 
in the hot plasma. 

We will consider some lattice gauge theory data 
concerning these peaks in the second half of this 
talk. 

3. Screening and J/i/; Suppression. 

One of the first theoretical suggestions for a 
signal of the existence of the quark gluon plasma 
concerned the heavy quark states J/'i/;. It was ar- 
gued and backed up with potential model calcu- 
lations, that quasi- free quark and gluon thermal 
screening would reduce the attractive forces be- 
tween the heavy quarks sufficiently to eliminate 
the binding energies and the J/-0 states them- 
selves In fact, the suppression of these states, 
again implied by the absence of the associated 
lepton pairs, is much more dramatic than that 
expected in hadron absorption models. 



6. Is there Nonperturbative, Semi- 
Classical Physics at Work in Heavy 
Ion Collisions? 

RHIC experiments suggest that there are early, 
strong, high multiplicity parton-parton collisions 
in heavy ion scattering processes. Perturbation 
theory, such as the QCD multiperipheral dia- 
grams 1^ which can fill the rapidity plot with a 
uniform density of gluons, fall far short of describ- 
ing the data. This suggests that larger length 
scales where the running coupling is in the inter- 
mediate coupling range are involved. This is the 
domain where semi-classical collective excitations 
of QCD are important. E. Shuryak has suggested 
that instantons and sphalerons px| ] may be caus- 
ing some of the surprising scattering events seen 
at RHIC. 
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7. What is a Lattice Theorist to do? 

We have taken a short survey into topics of 
interest to heavy ion physicists. Can lattice theo- 
rists have much impact here? Much of the physics 
here is stated in Minkowski space, natural to scat- 
tering processes. Employing Euclidean Lattice 
methods to tackle some of these questions will 
certainly take new initiatives and methods. As 
we will review in the second half of this talk, lat- 
tice methods do provide information on real time 
processes even within the present techniques we 
use. These include power spectra, Maximum En- 
tropy Methods and lepton pair rates. However, 
more comprehensive methods closer to first prin- 
ciples can probably be developed, too. 

Let's look at a few such topics. 

7.1. Lattice Approach to Nucleus Wave 
Functions 

The parton picture of QCD structure requires 
the use of the infinite momentum frame, or, 
equivalently light cone quantization. Each par- 
ton has a particular tranverse and longitudinal 
momentum in a wavefunction expansion, as writ- 
ten above. Infinite momentum frame Hamilto- 
nians have been written down in this language 
and two dimensional Hamiltonians which control 
the transverse dynamics have been studied. Lat- 
tice gauge theory versions of these systems can be 
constructed ||ll|] and equivalent three dimensional 
Lagrangian systems derived. 

One can hope that formulations of this sort 
could be used to address issues of parton satura- 
tion in the context of lattice methods which can 
handle honestly the intermediate coupling char- 
acter of these problems. The infinite momen- 
tum frame idea will eliminate vacuum fluctua- 
tions from the nonzero longitudinal fraction parts 
of the problem. 

In the context of continuum methods, there 
have been recent attempts along these lines by 
heavy ion physicists. An ambitious attempt to 
formulate deep inelastic scattering in Euclidean 
terms has been pioneered by O. Nachtmann jl3|. 
This looks like particularly fertile ground for a 
lattice incursion. 



7.2. Nonperturbative Physics for T > Tc 
and small /i 

As I will review in the second half of this talk, 
the region of the QCD phase diagram for small 
fi in the vicinity of the quark gluon transition Tc 
is accessible to conventional simulation methods. 
This is a recent development that opens up ex- 
perimentally relevant lattice calculations to the 
whole field. The thermodynamics of the forma- 
tion of the quark gluon plasma for T near Tc and 
small Baryon chemical potential /x should be stud- 
ied very professionally now. Some of the projects 
will be reviewed in the second half of this talk. 

We would certainly like to know how quickly 
the underlying physical mechanisms responsible 
for the transition from the hadronic phase to the 
quark gluon plasma phase change as /i is turned 
on for T w Tc. We know that there are substan- 
tial interactions in the quark gluon plasma just 
above Tc- How arc these mechanisms and how 
are bulk thermodynamic quantities effected as n 
is taken nonzero? 

The role of instantons in the formation and 
character of the quark gluon plasma needs more 
elucidation even at the /i = edge of the QCD 
phase diagram. According to the instanton liquid 
model ||l^, the instanton ensemble is disordered 
in the hadronic phase. This is essential for the 
model to produce chiral symmetry breaking and 
phenomcnologically reasonable numbers. It is be- 
lieved that as T passes through Tc, instantons 
and anti-instantons pair up and form "neutral" 
molecules. Chiral symmetry is restored in such 
an ordered ensemble. 

It would be interesting to verify this "binding- 
unbinding" picture of the phase transition by 
lattice methods and to provide estimates of the 
scales of the phenomena involved. This will not 
be easy because instantons and anti-instanton 
molecules will tend to "fall through" the lattice 
and be hard to identify, let alone measure. Im- 
proved actions will probably be needed. 

Such studies would be worth doing because we 
know from other simulations that there are sub- 
stantial interactions in the quark gluon plasma 
just above Tc- We would also like to know how 
instanton dynamics change as fi is turned on for 
T^Tc. 
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These are interesting questions which are ex- 
perimentaUy relevant. Conventional lattice meth- 
ods can address them now. 

7.3. Diquarks and Color Superconductiv- 
ity? 

Perturbative and instanton methods 
have been used to predict color superconductiv- 
ity for low temperatures and asymptotically large 
11. A large gap has been found as well as in- 
teresting symmetries (" color- flavor locking" p^). 
The large gap in the spectrum of the Cooper 
quark-quark pairs suggests that this state of mat- 
ter should occur in natural phenomena, such as 
the interiors of neutron stars, at high but not 
inaccessible densities. The lattice has been un- 
able to contribute to this field because of the sign 
problem in the Euclidean fermion determinant in 
SU{3) lattice gauge theory. 

The weak coupling, perturbative and instan- 
ton, approaches to color superconductivity have 
not addressed the numerical impact of confine- 
ment on their considerations, so estimates of a 
critical density are suspect. Much of the research 
has concentrated on SU{2) color which is quali- 
tatively different from QCD. Cooper pairs in the 
SU{2) model are color singlets and are degenerate 
with ordinary meson states because of the self- 
conjugate nature of quark and anti-quarks in this 
case. At chemical potentials greater than half the 
pion mass, the ground state of the model becomes 
a superfluid. 

Instanton enthusiasts have noted that their 
model provides attractive forces which favor the 
appearance of scalar 3 diquarks jl^]. The forces 
are only half as strong as the analogous ones in 
the SU{2) model where diquark states are degen- 
erate with mesons. This suggests one should find 
considerable scalar 3 diquark correlations in ordi- 
nary hadronic physics. To my knowledge, there 
are none. Perhaps lattice spectroscopists could 
scour their baryon wavefunctions and search for 
diquark correlations in their three body states. 

Aside from the issue of color superconductivity, 
it would be good to settle the old controversy con- 
cerning the possible viability of diquark models of 
baryon structure. 

Weak coupling descriptions of color supercon- 



ductivity don't need diquark correlation enhance- 
ments in ordinary hadronic states to be success- 
ful and accurate at nonzero chemical potentials. 
This is because the Cooper pairs of the BCS the- 
ory have very large, "floppy" wavefunctions which 
are much larger than protons ||l7[| . 

8. Contributions to Lattice 2002 

The parallel sessions at Lattice 2002 concern- 
ing QCD in extreme environments contained a 
number of interesting contributions. Because of 
the lack of space I will discuss just four topics 
that generated the most discussion. I will have 
to point the reader to the writeups of the parallel 
talks for details. 

The four topics will be : 1. The high T, low 
fi edge of the QCD phase diagram, 2. Thermal 
dilepton rates, 3. Real-time pion propagation in 
hot QCD, and 4. Four dimensional Four Fermi 
models at nonzero /i. 

8.1. High T, low /i edge of the QCD phase 
diagram 

Following the seminal work of Fodor and Katz 
last year the field has realized that the high 
r, low fj, edge of the QCD phase diagram is ac- 
cessible to a variety of mature lattice techniques. 
Fodor and collaborators presented arguments 
that their multi-parameter generalization of the 
Glasgow method has an efficiency which only falls 
off as a power, V^°' with a w 0.2 — 0.3, as the 
size of the lattice grows large, de Forcrand and 
Philipsen pO| have done simulations at imaginary 
chemical potentials using conventional Hybrid 
Molecular Dynamics algorithms and have found 
convincing, controlled results for fiq < ttTc/S. 
The Bielefeld-Swansea group have presented 
Taylor expansions in /i of the Fodor-Katz method. 
D.K. Sinclair |2^ simulated SU{3) lattice gauge 
theory at nonzero isospin chemical potential /x/ 
and found that Tdpi) resembles the curve Tdp) 
predicted by the other groups. 

Before we consider the first two contributions in 
more detail, recall some elementary points. The 
transition from the hadronic phase to the quark 
gluon plasma phase at T = Tc and /i = 0.0 
is actually a rapid crossover for the range of 
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quark masses presently under investigation. Most 
lattice theorists believe that it is also a rapid 
crossover for realistic quark masses, niu, rud and 
■nis, but larger scale simulations are needed here. 
As n is turned on one expects that the crossover 
Tc{n) to decrease slowly. However, for sufficiently 
large /i one expects the crossover to become a first 
order transition and it is this point " E" , the crit- 
ical endpoint, Te and fj,E, that one would like to 
know [ p^ . Present simulations indicate that /i^ 
is quite large However, those simulations 

are on small lattices with relatively large quark 
masses. They should be repeated on large lat- 
tices with realistic quark masses to get a relevant 
answer for Te and fj,E- if f^E turns out to be 
very small, there would be dramatic effects, po- 
tentially, in the projectile fragmentation regions 
at RHIC. 
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Figure 2. Evolution of the crossover to the end- 
point of a first order line of transitions as hb in- 
creases. 



Fodor and Katz [Q introduced multi- 
parameter reweighting to improve the overlap of 



lattice configurations on the edge of the phase 
diagram with particular points inside the phase 
diagram, 

Z(m, fi,(3) ^ / D[/e-'^''°=('^"'^^ det M(mo, /i = 0, C/) 



det M (m, fj,, U) 
detM(mo,/i = 0, U) 



} 



This formula organizes the calculation of the 
Partition function at a point inside the phase 
diagram, by making configurations on the edge 
of the phase diagram (first line) and reweight- 
ing (second line) in two variables, /i and (3. The 
idea is that configurations at Tc experience both 
the hadron and the quark gluon plasma phase 
so there should be a considerable overlap with 
configurations nearby but at different /i and (3 
values. Care must be taken here : as /i is in- 
creased, (i must be increased slightly so the sys- 
tem remains near the crossover coupling for that 
/i value, /9c(m)- As ^ and (3 are increased, one 
monitors the acceptance rate and keeps it near 
50 percent . In this way lines essentially par- 
allel to /?c(m) but displaced into the hadronic or 
the plasma phase have also been investigated and 
the equation of state has been calculated for a 
range of T and /z values. They have shown that 
the crossover sharpens as /i grows and the end- 
point "E" is approached The position of the 
endpoint was found in last year's contribution but 
studying the finite size scaling properties of the 
Lee- Yang zeros of the Partition function p8[ |. 

Following the paths in phase space where the 
authors monitor the efficiency of their overlap 
method, they are now computing the Equation of 
State . The crossover is predicted to sharpen 
considerably as fiB increases. 

The authors argue that at least for volumes up 
to 12^ X 4, they can simulate chemical potentials 
up to /I oc V' with 7 w 0.2 — 0.3. And they 
argue that the efficiency of their multi-parameter 
overlap method scales as V~°' with a w 0.2 — 0.3, 
so the method should allow fairly large lattices 
to be simulated before its inefficiencies become 
overwhelming. 

The reader should consult the parallel session 
contribution by these authors for plots and details 
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, new method 



Glasgow method > 



1 critical hne 




Figure 3. Paths in the phase diagram used by 
Fodor and Katz compared to Glasgow's original 
method. 



|19|. 

The second contribution presented on this topic 
by de Forcrand and Philipsen simulated lat- 
tice at imaginary chemical potential using the 
conventional hybrid molecular dynamics algo- 
rithm. It followed the crossover, /3c(m)> by com- 
puting the plaquette susceptibility and monitor- 
ing its peak as a function of f3 for fixed values of 
the imaginary chemical potetial. 

The authors argue that the curve /3c (m) is an 
even, analytic function of /i. Writing /3c (/i) — 
^,^c„(a^)^", they find that only the first two 
terms are needed to describe the data accurately. 
The analytic continuation to real fj, is, therefore, 
trivial. Their resulting curve is in agreement with 
that of Fodor and Katz, although they haven't in- 
vestigated the order of the transition as a function 
of fi yet. 

The imaginary chemical potential method is 
very conservative and effective. It will work on 
large lattices with the efficiency of an ordinary 
hybrid molecular dynamics code. It can only 



treat small chemical potentials because imaginary 
chemical potentials induce tunnelling to other 
2' (3) vacua when the imaginary chemical poten- 
tial equals TrT/3 ||]. For T near Tc w 170 MeV., 
the method is limited to /i < 500 MeV. 

The reader is referred to the relevant contribu- 
tion at the parallel session for quantitative details. 
In addition, similar work has appeared by D'Elia 
and Lombardo p^ . 

We can look forward to considerable progress 
in understanding this strip of the QCD phase di- 
agram now that several simulation methods and 
groups are involved. 

8.2. Thermal dilepton rates 

Dilepton rates were discussed above as an in- 
formative probe into the dynamics of the quark 
gluon plasma. This topic also provides a nice 
example of traditional lattice gauge theory con- 
fronting new experimental scattering data. 

The dilepton rate, which is a function of fre- 
quency uj and momentum p, can be written in 
the form lEl, 



dW 



1 



-(Ty(w,_p;T) 



(3) 



dwdp 277r2 w2(e"/^ - 1) 

where the spectral function av{^iP]T) is re- 
lated to the current-current correlation function, 
Gy(T,p;T) = /dxe^P - < J(;(t, f)4^(0, 0) >, 
on the periodic lattice as. 



r ( -rj.. /■ , I ^^, cosh(u;(T-l/2T) 
J smh([x)/2J) 

On a X Nr lattice, one can measure the cor- 
relation function Gy at each "time" slice, t^T = 
k/Nr, = 1, 2, Nt and then attempt to invert 
the relation between the spectral function and the 
correlation function. 

There are two approaches to executing this in- 
version. First, one might use a model dependent 
functional form which depends on several param- 
eters and determine (T\/(w,p;T). This approach 
has had only limited success. Alternatively, one 
might use the Maximum Entropy Method (MEM) 
p^ ] and find the "most probable" spectral func- 
tion. Since the lattice simulation gives the cor- 
relation function only at a discrete set of points 



(4) 
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Tfe, each with statistical errors, the problem of 
finding the continuous function cry of the contin- 
uous variable lu is ill-posed. The Maximum En- 
tropy Method can only produce a "most proba- 
ble" answer with a measure of its reliability. The 
method, however, is well-known and quite suc- 
cessful in optics, where it is used for image recon- 
struction. Again, the reader should consult the 
parallel sessions for a wealth of examples that il- 
lustrate the method's strengths and weaknesses. 

At past conferences the Bielefeld group has pre- 
sented data 1 25 1 which have shown that the p and 
4> peaks in the dilepton spectra are suppressed in 
the presence of the quark gluon plasma, in quali- 
tative accord with the experimental trends. How- 
ever, the dilepton signal is still somewhat larger 
than that found for free but hot quarks and glu- 
ons, at least for lo/T between 4 and 8. These 
simulations were done above the transition, at 
1.5Tc and 3Tc, so measureable attraction in these 
light quark states at these high temperatures is 
somewhat surprising. Results such as these sug- 
gest that although bulk thermodynamic behav- 
ior of the quark gluon plasma might be close to 
the Stephan-Boltzmann limit, detailed hadronic 
structures are not completely washed out. 

At this conference S. Datta ||2^ presented 
Bielefeld simulation results at 0.9Tc and 1.2Tc for 
the charmonium states. The S-wave states had a 
diminished, perhaps by a factor of 2, signal above 
the transition where potential models predicted 
there would be no signal at all. The P-wave sig- 
nal, however, was reduced by a factor of about 
7 across the transition suggesting that screening 
has removed these charmonium states from the 
spectrum. 

8.3. Real-time pion propagation in hot 
QCD 

A study of Effective Lagrangians by Stephanov 
and Son |Q has uncovered curious features in the 
propagation of pions as the quark gluon plasma is 
approached by heating hadronic matter, T < Tc- 
In fact, the pion dispersion relation is found to be 
determined by temperature dependent features of 
the system which are accessible to conventional 
Euclidean lattice simulations. The three features 
one needs are 1. the pion screening mass, 2. the 



pion decay constant, and 3. the axial isospin sus- 
ceptibility. 

This work constitutes a nice, informative ex- 
ample where "real-time" quantities are accessible 
to static simulations. Of course, the most popu- 
lar example of this is the relation of the speed of 
sound Vs in a system to its bulk thermodynam- 
ics, v1 = dP/de, where P is the bulk pressure 
and e is the energy density. Other examples ex- 
ist in condensed matter physics, as well. For ex- 
ample, the propagation features of spin waves in 
anti-ferromagnets are determined in large part by 
static correlation functions. 

These authors determine the real part of 
the soft pion dispersion relation. 



u {p -f m ) 



(5) 



where we identify u, the pion's speed, to, the 
pion's screening mass, and Wp=o = um, the pion's 
"pole mass". The authors derive three low en- 
ergy theorems at nonzero T. First, v? = /^/x/5, 
where / is the pion decay constant and x/5 is 
the axial isospin susceptibility. These last two 
quantities are related by the other two theorems, 
which represent the generalizations of the Gell- 
Mann Oakes Renner relation to finite tempera- 
tures, pw? = = —f^q < i)ip >■ Since 
these quantities are accessible to conventional lat- 
tice simulations, we have a nice illustration of 
determining "real time" physical quantities from 
static, Euclidean simulations. 

The strategy of this development is interesting 
and quite general. One begins with the micro- 
scopic quark Lagrangian underlying lattice QCD, 



T3 



(6) 



as well as the low energy effective Lagrangian that 
describes the properties of the low energy mesons. 



f2 f2 

L^ff = if TrVoSVoSt - ij-TrdiY.diT.'^ 



1 



-m^f^ReTrY. 



(7) 



where E = exp(iT"7r''/F^), and VqS = doT. - 
^/i/5(T3E + Sra). The Effective Lagrangian 
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should be familiar from strong interaction phe- 
nomenology and its extension to nonzero chemical 
potential is determined by the symmetries of the 
underlying microscopic Lagrangian without any 
numerical or phenomenological ambiguity [29|] . 

By matching the predictions of the microscopic 
theory with those of the low energy Lagrangian, 
the authors determine their low energy theorems. 

Note how curious pion propagation is as T ap- 
proaches Tc- The pion propagation speed ap- 
praoches zero, as does the pion decay constant 
/. A further investigation of the scaling laws gov- 
erning the approach to the quark gluon transition 
reveals that u vanishes faster as T approaches Tc 
than the screening mass m diverges so the 
pole mass nip = um approaches zero as T ap- 
proaches Tc- Since the pion's speed vanishes in 
this limit while its pole mass does also, perhaps 
we can call the pion an "infinitely massive mass- 
less" state! Hopefully, some of these features will 
show up in final state studies at RHIC. 

8.4. Four dimensional Four Fermi models 
at nonzero fi 

As advertised in the first part of this talk, it 
would be good for all concerned if there were 
more interaction between the lattice community 
and the phenomenologists interested in the quark 
gluon plasma and in color superconductivity. One 
step in this direction has been taken by the 
Swansea group ]30| which has found diquark con- 
densation in strongly cutoff 3 + 1 dimensional 
Nambu Jona-Lasinio models, of the sort discussed 
by Nuclear theorists, instanton phenomenologists 
and others. The value in the lattice simulations 
is, as usual, that their results are exact in princi- 
ple, because they include all the fluctuations. 

The Lagrangian which has been simulated is 
just one of the original Nambu Jona-Lasinio mod- 
els, 

J{i'^CJ,T2e^,i;f)] (8) 

where the last two terms are source terms for 



diquark condensation which are taken to zero 
U ~^ 0) at the end of a simulation meant to cal- 
culate the diquark condensate. These terms play 
the same role as a small bare fermion mass term 
in a lattice simulation of the chiral condensate — 
the simulation must pick out a preferred direction 
in the space of the condensate, calculate the con- 
densate and see if it survives in the limit where 
the source is taken away. In fact, simulations of 
this model indicate that as ^ increases, a diquark 
condensate appears in the system and its fi depen- 
dence and magnitude are close to the predictions 
of the leading term in the large N expansion. The 
reader should consult the parallel sessions of these 
conference proceedings for plots [ pO[ . 

One notable feature of this calculation is that 
it is the first lattice simulation in 3 + 1 dimen- 
sions with a traditional fermi surface at nonzero 
chemical potential. Future work in this direction 
will be the study of additional four fermi mod- 
els suggested by the instanton liquid model. One 
complication will be the cutoff nature of these 
models — these models are logarithmically trivial 
and they are useful only when strongly cutoff, 
with their parameters adjusted so that observ- 
ables such as the pion decay constant and the 
chiral condensate take on physical values. Luck- 
ily, the Nuclear physics community has used such 
models for some time and know the "best" way 
to deal with these issues. The lattice model may 
reach an impasse, however, because when com- 
plicated four fermi interactions are entertained, 
it becomes impossible to avoid complex fermion 
determinants. In addition, these models do not 
confine quarks in their low temperature phase, so 
they are subject to criticisms discussed earlier. 
Nonetheless, these studies should be informative 
and represent a test of some of the analytic, ap- 
proximate approaches to the high /i, low T super- 
conductivity phase. 

9. Conclusions 

Lattice gauge theory has had much to say 
about the finite temperature behavior of QCD 
because it can handle the three length scales of 
QCD, perturbative behavior at weak coupling 
and short distances, semi-classical coherent be- 
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havior at intermediate couplings and intermedi- 
ate distances, and confinement at strong cou- 
plings and large distances. It is now able to study 
the variable T-low strip of the phase diagram 
and make predictions relevant to heavy ion ex- 
periments. I hope that its lessons will lead to 
additional formulations and strategies to study 
real time physics relevant to collisions and nu- 
clear structure in the not-too-distant future. 
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